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Abstract 

The Schrodinger picture description of vacuum states in curved space- 
time is further developed. General solutions for the vacuum wave func- 
tional are given for both static and dynamic (Bianchi type I) spacetimes 
and for conformally static spacetimes of Robertson- Walker type. The for- 
malism is illustrated for simple cosmological models with time-dependent 
metrics and the phenomenon of particle creation is related to a special 
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form of the kernel in the vacuum wave functional. 
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1 Introduction 



Many of the unusual and potentially paradoxical phenomena which occur in quan- 
tum field theories in curved spacetimes reflect the nature of the vacuum state, for 
a review see for example P], EL EL M. However, whereas in Minkowski spacetime 
we have a simple and unambiguous prescription for determining the vacuum, the 
specification of the vacuum state in a general curved spacetime involves many 
subtleties. Indeed, it is not assured that a state satisfying the defining attributes 
of the Minkowski vacuum even exists for a general spacetime. 

In this and a companion paper ||, we discuss the nature of the vacuum 
state for a broad class of flat and curved spacetimes and develop a formalism 
of sufficient generality to allow a future application to the vacuum structure of 
black hole spacetimes and the associated Hawking radiation || |J . 

This formalism is based on the Schrodinger wave functional. The essential 
reason for this choice is to avoid the conventional description of the vacuum 
as the 'no-particle' state in a Fock space of states generated by the creation 
operators of particles defined with respect to a particular mode decomposition of 
the quantum field. This approach is inherently problematic in curved spacetime, 
where there is no uniquely favoured mode decomposition and no guarantee that 
the usual concept of a particle is a good description of the spectrum of the theory. 
The Schrodinger wave functional circumvents this second problem by giving an 
intrinsic description of the vacuum without reference to the spectrum of excited 
states. 

The wave functional that defines the vacuum state satisfies a functional Schro- 
dinger equation, which presupposes a foliation of spacetime into a series of space- 
like hypersurfaces, the Schrodinger equation describing the evolution of the wave 
functional between successive hypersurfaces. The wave functional therefore de- 
pends on the choice of foliation. | Our expectation is that physical quantities 
should be independent of this foliation. However, intermediate quantities such 
as Green functions may have a foliation dependence. This dependence should 
be controlled by a functional Ward identity (somewhat analogous to the depen- 
dence of Green functions on the choice of gauge in ordinary quantum field theory). 
Physical quantities should depend only on the spacetime manifold and boundary 
conditions. 

3 Notice that the choice of foliation is quite distinct from a choice of coordinates. The wave 
functional is independent of any particular coordinates used to describe either the spacetime or 
the foliation (see section 3). In the Schrodinger formalism, the choice of foliation replaces the 
mode decomposition in the conventional approach. One of the advantages of this formalism is 
the clear separation it provides between the choice of vacuum state and the issue of coordinate 
independence. 
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In this paper, which is intended as the first step in a general programme, we 
consider a class of globally hyperbolic spacetimes, viz. manifolds with the topol- 
ogy Ai = 1Z x E where E are spacelike hypersurfaces. Within this category, we 
distinguish manifolds which are 'static' and 'dynamic', essentially implying that 
the metric depends only on the spacelike and timelike coordinates respectively 
The class referred to as 'dynamic' is in fact the Bianchi type I spacetimes. We 
also consider more general conformally static spacetimes where the conformal 
factor is purely time dependent. This includes the cosmologically important class 
of Robertson- Walker spacetimes. 

The ultimate goal of our work is a full analysis of the physics of the vac- 
uum state in the Schwarzschild and Kruskal spacetimes describing collapsed and 
eternal black holes. In this and the companion paper, we develop all the nec- 
essary theoretical formalism and, as a warm-up, illustrate the techniques on a 
number of simple illustrative examples. The solutions for static and dynamic 
spacetimes are sufficient to encompass the full Kruskal metric and we check ex- 
plicitly that subtleties involving foliation dependence and the interpretation of 
vacuum states are under control. The extension to black hole spacetimes is in 
principle a straightforward application of these techniques, although the presence 
of real event horizons separating the static and dynamic regions is an important 
new feature. 

A number of the results and techniques described here can be found already 
elsewhere in the literature. In these cases, we have referenced the original works 
as fully as possible, but have included the material as we have been concerned to 
present a coherent and systematic study of the Schrodinger wave functional. 

The paper is organised as follows. In section 2 we review the Schrodinger 
picture formulation of quantum field theory in Minkowski spacetime and derive 
the Schrodinger equation and its solutions for the vacuum and excited states for 
a free scalar field. We also discuss some properties of a general Gaussian state 
which will be needed in sections 4 and 5. 

Then in section 3 we construct the Schrodinger wave functional equation for an 
arbitrary foliation of a globally hyperbolic spacetime with topology M. = 1Z x S. 
This makes explicit the dependence of the equation, and therefore the vacuum 
wave functional, on the 'lapse' and 'shift' functions specifying the foliation. The 
simplification for special coordinate choices which mirror the foliation is shown. 
We find explicit solutions for the vacuum wave functional for the special cases 
of static and dynamic spacetimes and a particular class of conformally static 
metrics including the Robert son- Walker spacetimes. This analysis shows clearly 
that while there is frequently an essentially unique vacuum state for a static 
spacetime, in the dynamic case there is necessarily a one-parameter family of 
vacuum solutions of the Schrodinger equation. 
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In sections 4 and 5, these techniques are illustrated using two simple two- 
dimensional 'cosmological models' which have been much discussed in the liter- 
ature. The first represents an asymptotically static universe that undergoes a 
period of expansion while the second describes a universe which expands from an 
initial singularity before recollapsing. These are of particular interest as models 
illustrating cosmological particle creation. We solve the Schrodinger wave func- 
tional equation exactly to determine the candidate vacuum states and discuss 
carefully the interpretation that particles are created by the expansion. 

A number of technical results on solutions and Green functions for various 
differential equations occurring in the text are given in the appendices. 



2 Schrodinger Picture in Minkowski Spacetime 

To introduce the Schrodinger formalism, we begin by reviewing the wave func- 
tional equation and its solutions in Minkowski spacetime. (For a review see, e.g. 
[0.) Consider a massive scalar field ip(x) with action^ 



S = I d d+1 x C (1) 

where the Lagrangian density is 

£ = ^{^(^)(M-™V} (2) 

The canonical momentum conjugate to the field is defined by 

dC 

vr(x) = — = WX) (3) 

The Hamiltonian is constructed by a Legendre transform of the Lagrangian 

H(tt, if) = J d d x [7r0 - C] (4) 

= \ Jd d x {tt 2 - + m\ 2 ) (5) 

The system can now be canonically quantised by treating the fields as oper- 
ators and imposing appropriate commutation relations. This involves the choice 
of a foliation of the spacetime into a succession of spacelike hypersurfaces. In 
Minkowski spacetime, it is natural to choose these to be the hypersurfaces of 
fixed t and impose the equal-time commutation relations 

[(p(x,t),7r(y,t)} =i5 d (x-y) 



4 Our conventions are as follows. For generality we consider a spacetime of d+ 1 dimensions. 
The metric signature is (+,—,..., —). Coordinates for the spacetime point x are denoted either 
by (a; , x l ), where i = 1, . . . , d, or by (t, x). 
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[<p{x,t),ip{y,t)\ = [n(x,t),n(y,t)} = (6) 

In the Schrodinger picture, we take the basis vectors of the state vector space 
to be the eigenstates of the field operator ip(x) on a fixed t hypersurface, with 
eigenvalues <f)(x), 

<p(x)\(/>(x),t) = <t>(x)\(/>(x),t) (7) 

Notice that the set of field eigenvalues <p(x) is independent of the value of t la- 
belling the hypersurfaces. This will be especially important in the curved space- 
time context. 

In this picture, the quantum states are explicit functions of time and are 
represented by wave functionals ^[<j)(x),t]. Operators 0(ir,(p) acting on these 
states may be represented by 

o(*(x)M*)) ~ °(- l J^uy^) ^ 

consistent with the commutation relations (||). The Schrodinger equation, which 
governs the evolution of the wave functional between the spacelike hypersurfaces, 
is 




+ m 2 2 }*[<M] (9) 



To solve this equation, we make the ansatz that, up to a time-dependent 
phase, the vacuum wave functional is a simple Gaussian. This is in direct analogy 
with the standard solution of the harmonic oscillator in quantum mechanics. We 
therefore write 

V [<f>,t] = N (t) exp{-iy>x J d*y <!>(£) Gfey)^} = N (t)M<P] (10) 
Substituting into (^) and comparing coefficients of [0°] determines N (t) from 

^ = 4/4GM (ii) 

while comparing coefficients of [4> 2 ] gives an equation for the kernel G(x,y): 

d d z G(x, z)G(z, y) = (- V 2 + m 2 )5 d (x - y) (12) 

The kernel equation can be solved by using several methods, either by working 
in momentum space, or coordinate space, or using the Schwinger-DeWitt repre- 
sentation of the propagator ||, ||. 
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We show here the momentum space approach, which is the most convenient 
in Minkowski spacetime. Introducing the Fourier transform of the kernel, 

G(x,y)= J 0^- d e G(k) (13) 
the kernel eq. (|12D reduces to 

G(k) 2 = k 2 + m 2 (14) 
The inverse of the kernel, A(x,y), is defined as 

J d d z G(x,z)A{z,y) = 5 d (x-y) (15) 
and is therefore given explicitly by 

A/ , r d d k j^yJ 

Afe ^/(SF^ (16) 

The kernel is simply related to the Green functions of the theory (see appendix 
B). A general Green function can be written as 

f d d+1 k e'- ) ~ ik (x -yo) 

g(x > y) = J Wy^ — — (17) 

where u 2 = k 2 + m 2 . The ko integral has poles at ko = ± uo and can be evaluated 
with different boundary conditions, determining how the contour goes round the 
poles. In particular, the Wightman function Q + (x,y) is defined by choosing the 
contour integral to close around the k = uj pole as in Fig(|3]) of appendix 0. For 
equal times, this Green function reduces to 

„ , 1 f d d k e^-v) t 

G + {x,y)\ E T = - / , (18) 

2 J ( 27r ) ^Jk 2 + m 2 
_ j_ 

2 

Kd_i(m\x-y\) (19) 



m 



2ir \27r|x — y\J 

The inverse kernel is therefore simply twice the Wightman function evaluated at 
equal times, i.e. 

A(x,y) = 2g + (x,y)\ ET (20) 

The reason for this identification is made clear in appendix where we discuss 
the representation of vacuum expectation values and two-point functions in the 
Schrodinger picture. 
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The kernel itself is 



G(x,y) = \g+(x,y)\~Bh- (21) 
= /705V / i r ^^ , ^ ) (22) 



d , 1 
2 + 2 



= - 2 (ss=iJ w^-tf < 23 > 

The time dependent phase factor in the vacuum wave functional is 

N Q (t) = e~ iEot (24) 

where 

E = 2/ ( 25 ) 

= 1 y \/^ 2 + ™ 2 <^(°) ( 26 ) 
is the (divergent) ground state energy. 

There are several ways to derive the excited state solutions to the Schrodinger 
wave functional equation, but with a view to the curved spacetime generalisation, 
we start with a method that does not presuppose a particle interpretation of the 
excited states. 

Making the following ansatz for the first excited state 

* 1 [4>,t]=2N 1 (t) [d d x[d d y(f>(x)G(x,y)f(y)M<i>} = (27) 



the Schrodinger equation is solved with the same kernel (|i2[ ) as for the ground 
state, provided 

Nt(t) = e~ iElt (28) 
where the energy gap AEi — Ex — E satisfies 

AE X J d d y G(x,y)f(y) = (-V 2 + m 2 )f(x) (29) 

This is readily solved if f(x) is an eigenfunction of the operator (—V 2 + m 2 ). We 
therefore choose fk(x) = e 1 --, which has an eigenvalue of (fc 2 + m 2 ). We then 
find 

A Ei [ d d y G{x,y)f k {y) = (k 2 + m 2 )f k (x) 



which reduces to 

AEt G(k) = k 2 + m 2 
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so the energy gap is 



AE 1 



k 2 + m 2 



UJ 



(30) 



Clearly this admits a particle interpretation where the first excited state is 
a one-particle state and E\ is its energy. This is apparent from the energy gap, 
AEi = \Jk? + m 2 , which is the energy of one particle with momentum k. The 
wave functional carries a label k, introduced through the eigenfunction f^, which 
specifies the momentum of the particle. 

The Schrodinger equation can also be solved directly (as in the quantum 
mechanical harmonic oscillator) for an n-particle state, giving 



(31) 



where the H n are generalised Hermite polynomials defined using the kernel and 
eigenfunctions fk. contracted with the field <fi as in fl2~7|) . The energy of the nth 
excited state is simply 

n 

E n = E + J2^(ki) (32) 
i=i 

where k^ is the momentum of the ith particle. 

The excited states can also be described using creation and annihilation op- 
erators. We define 



a(k) 



d d x 
d d x 



d d y G{x,y)(f>{y) + 
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6<j>(x) 



UJ(j)(x) + 



6 



ft{x) 



and 



a f (fc) 



d d x 
d d x 



Jd d yG(x,y)^y) 
5 



5(f)(x) 



fk(x) 



UJ0[X I — 



S(f)(x) 



consistent with the commutation relations 

[a(k) , a j (p)} = 2 (2ir) d to 5 d (k - p) 



(33) 
(34) 

(35) 
(36) 

(37) 



These can be used iteratively to construct all the excited states, starting from 



d(k)ip = 



^(k)^ = if)! 



The nth excited state is labelled by the n particle momenta, viz. 

ipnih, ...,k n )= a^) . . . a j (k n )ip 



(3? 



(39) 
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Orthogonality relations between these states follow directly from the commu- 
tation relations of the creation and annihilation operators. We have 



(Mkl,-..,k n Mm(P 1 ,---,pJ) = <U{ n2^(27r) d 5%-£.) 

i=l 

+ ...n\ permsj^ol^o) (40) 

where the inner product is defined by 

(tfxltfa) = [V<P *5(0)* 2 (0) (41) 



and uJi = sjki + m 2 . The Hamiltonian and particle number operators are written 
in terms of a and a) as usual as 



d d k 



and 



A J (2tt) £ 



AT 



d d fc 1 



(42) 



(43) 



(2vr) d 2cj 

It will often be convenient to express the wave functional for the vacuum 
and excited states in terms of the Fourier transform <f)(k) of the field eigenvalues 
defined in eq.(0). We define 

and since 4>(x) is real, <f)*(k) = <p(—k). The vacuum wave functional is therefore 

1 r d d k 



exp 



(45) 



2 J (2vr) d j 
In terms of <f)(k), the creation and annihilation operators simplify: 

a(k) = ujm + -J— (46) 
d(p* (k)_ 

a\k) = uo0*(k)-JL- (47) 



5<j)(k)_ 

where 5(j)(p) / 5(j)(k) = (2ir) d 5 d (k - p) and 50*(p_)/50(A;) = (2tt)^(A; + p) . The 
excited states for particles with momenta k { are therefore given by 



2[u; 1 0*(& 1 )]V>o 



(4* 



V> 2 = {4[ Wl 0*(i 1 )][^ 2 0*(A; 2 )]-2cu 1 (27r) d (5 d (A; 1 + A; 2 )}^o 



(49) 
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etc. The general n-particle wave functional is expressed in terms of generalised 
Hermite polynomials with n momentum labels, viz. 



-iE n t 



^n[</>;ii,---,in]^o[0] 



[f] 

= e~ iEnt J2 {(-l) m 2 n " 



m=0 



n^i(27T) d 5(A; 2 ,_ 1 +^) 



.i=i 



n ^*(^) 

j=2m+l 



+ ...(2m-l)!!( 2 ^)perms}^ (50) 



where there are (2m — 1)"Q^) permutations of the pairing of momenta at each 
stage of the sum. The sum is from zero to | - which is the largest integer < |. 
Products Ili=o with b < a are defined to be 1. 

All this analysis is based on the specific foliation of Minkowski spacetime in 
which the family of spacelike hypersurfaces £ are equal-time hypersurfaces and 
the Hamiltonian evolution is along the timelike Killing vector field t|. While this 
is convenient, it is by no means a unique choice. In section 3, we set up the 
Schrodinger equation for an arbitrary foliation in curved spacetime. This formal- 
ism may of course be specialised to Minkowski spacetime by simply replacing the 
metric by r]^ throughout. 



2.1 General Gaussian state 

In the discussion of the cosmological models in sections 4 and 5 we encounter 
Gaussian states where the kernel is time dependent and not simply equal to u, 
the eigenvalue of the wave operator. In anticipation of this, we show here how 
such a state in Minkowski spacetime may be expressed as a superposition of 
excited many-particle states. 

Consider a solution of the Schrodinger equation of the form 

*ftf] = JV(f)expj~/ 0- d m,t) |0(£)| 2 } (51) 

where 

N(t) = exp 
and Q(k; t) is a solution of 

.dfl(k;t) 2 2 

% — uj — \l 

dt 
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This can be written as a linear superposition of n-particle states as follows: 

d k/i 



n=0 



n 

i=l 



(27T) 



c n(^Il; • • • ih.n)' l Pn\h.l 



i • • • ) iiin ) 



(54) 



where the (time-dependent) expansion coefficients are 



c n (hi 



, . . . , n ) 



, . . . , fc n )^ 



(55) 



and we are assuming constant normalisation factors are included such that (^ol^o) 
and (v&lvP) are both 1. (For the justification that this is possible for the general 
Gaussian state, see section 3.) The n-particle states have been given earlier. We 
therefore have 



V^niti, ■ ■ ■ ,kn)* = I V$H n foh, . . . ,k n ] e-^^ {n+u3) im{2 (56) 



The odd moments of the integral vanish and therefore the coefficients c n vanish 
for n odd. For even n, we find 



Cn{h\i ■ ■ ■ i kn) — \~ 

1 n 



n/2 

n 



27T) d 5 d (k 2i _ 1 + k 2l ) /W-i-fi 



2i-l 



+ . . . {n — 1)!! perms . . . j^ol^} 



(57) 



where f2« = Q{ki'i t) and the (n — 1)!! permutations specify different combinations 
of the paired momenta. 

It is particularly interesting to look at the expectation value in the general 
Gaussian state of the number operator N(k), which counts the number of particles 
with prescribed momentum in the expansion (E3). This is 



[®\N{k)\^) 



1 

1 f 



[V\tf{k)a(k)\V) 



- (2n) d 5 d (0) 



(58) 



The expectation values can be easily evaluated using techniques described in 
appendix C. For complex Q, we find the expectation value of the number density 
to be 

(u -Q)(uj -Q*) 



(59) 



This result can also be recovered using the explicit expressions for the expansion 
coefficients c n given above. 
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3 Schrodinger Picture in Curved Spacetime 



We consider globally hyperbolic spacetimes of topology M. = 1Z x E with a 
global timelike vector field. To provide a general foliation of the spacetime, we 
use the formalism of embedding variables |TD], [IT], [Tj|. We therefore choose a 
family of spacelike hypersurfaces E, with intrinsic coordinates £ l , labelled by a 
'time' parameter s and consider evolution along the integral curves of J^. The 
embeddings characterising the foliation are maps x : S — > Ai which take a point 
on the hypersurface E to a spacetime point x M = 2 M (s,£ l ). 

First we need to construct the projections of spacetime quantities normal and 
tangential to the hypersurface E. The tangential projections are defined using 

The normal to the surface, n M , is timelike and is uniquely defined by 

= (61) 

with g^n^rty = 1. Any spacetime tensor can be expanded in terms of its projec- 
tions normal and tangential to E. We can also introduce an induced metric on 

hij = 9^PiP U j (62) 

and raise and lower hypersurface indices, e.g. define p 1 ^ = h^p^g^. 
Now consider the deformation vector of the foliation, defined as 

os 

= Nn^ + Wp? (63) 

where N and N l are the 'lapse' and 'shift' functions of the embeddings given by 

N = ri^N^ 1 (64) 
N* = p^JV" (65) 

In terms of these functions, the metric interval can be written as 

g^dx^dx u = (N 2 + WN^ds 2 + 2NidsdC + hijdCd^ (66) 

We also need the inverse relations 

= N (6?) 

•j... AT* 

+ Pa (68) 



~dxv ~ TV 
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and the Jacobian 



det 



dx^ 



(69) 



With these preliminaries complete, we can now construct the Hamiltonian 
describing the evolution appropriate to this foliation, i.e. along the integral curves 
of the timelike vector field J^. We consider a massive scalar field theory coupled 
to gravity. The action is 

s = -f d™xv=g\<rp-p- 

2Jm l dx^dx u 

and so in terms of the embedding variables, 



(m 2 + iR)^ 



(70) 



S 



- [ ds [ d d £det\J^- 
2Jr h - W(s,¥ 



dx^ dx v 



+2 



ds dC . 



where tp = |^ and <9jy? = J^. Using eqs.(|6T|), ( |68| ) and (|69"D, the action can be 
written in the form 

S = J R ds J^d\ £ (72) 
where the Lagrangian density is, 



+ h 13 



(d lV )(d 3 p)-(m 2 + ^R)p 2 } (73) 



The conjugate momentum is defined as usual to be 



dip N 

The Hamiltonian is given by the Legendre transform 

H = Jj\[ap-C} 



f d d i INK + WUi 



where 



H = Iv^{-^-^(^)(^) + (m 2 + ^V 2 } 



n 



(74) 

(75) 
(76) 

(77) 
(78) 



Canonical quantisation is implemented by imposing equal-s commutation re- 
lations for the field operators on the hypersurface E, viz. 



i^,s),n((,s)] = i 6 d (Z-0 
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[</?(£, s), <p({, s)} = [tt(£, s), tt(C, s)] = (79) 

In the Schrodinger picture, the quantum states are represented by wave function- 
als ^f[</>(^),s;N,N % ,hij] where the variables 0(0 are the eigenvalues of the field 
operator on the equal-s hypersurfaces, i.e. 

<p(x)\<j>(0,s) = 0(010(0, s) (80) 

As before, the commutation relations are realised by representing operators in 
terms of the field eigenvalues 0(0 and functional derivatives viz.f] 



The Schrodinger wave functional equation gives the evolution of \I/[0(O, s] : 




This equation holds for an arbitrary foliation (specified by N, N l and hy) and 
makes explicit the foliation dependence of the wave functional. In practical ap- 
plications, however, it is often convenient to choose spacetime coordinates which 
reflect the desired foliation. In this case, we simply identify the embedding vari- 
ables (s, with the spacetime coordinates (t,x). The lapse and shift functions 
are then N = y^oo an d N l = (so that goi = 0) while the induced metric is just 
hij = Qij. The Schrodinger equation reduces to 

l ~dt = 2J dd - v ^ 9 VJW ~ gij{di(t)){dj4,) + {m2 + ^ 2 }* (83) 

This is the form we shall usually use in the remainder of this paper. 

We can search for solutions to this Schrodinger equation in analogy to those 
found in Minkowski space. For the vacuum state, we make the ansatz 

^o[<P(x),t}=N (t)M<P,t} (84) 

where 

MM = ex p{~2 / d< W~^ j &v\pK t&G&iDWty)} ( 85 ) 

This has the same form as before, except that the kernel may now depend on the 
'time' t since the metric is in general time dependent. The spatial integrals in- 
clude the correct measure to ensure that they are invariant d dimensional volume 
elements. 

5 The functional derivative is defined so that 6<p(y) / '8<fi(x) — 5 d (x — y). The delta function 
density is given by S d (x,y) = W—h x )~ 1 5 d {x — y) and satisfies J d d x_y/—h x 5 d (x, y)f(x) = /(?/)• 
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Substituting this ansatz into the functional Schrodinger equation results 
in the time dependence equation 



dlnJVo(t) 
dt 



l -J d d x y/-h xy /gfo G(x, x; t) 



m) 



and the kernel equation 




i-^ {y/hxhy G(x,y;,tf) = J d d z \[^h z \[g^ \jh x h y G(x,z;t)G(z,y;,t) 



- Jh x h y Jg§ (Pi + m 2 + £R) X 5 d (x, y) 



(87) 



where = -^^di(g 1 ^ —gdj) is the spatial part of the Laplacian (see appendix 
A). Notice that this assumes that on the boundaries of the hypersurface £ the 
fields or their derivatives vanish, since the surface terms in an integration by parts 
have been omitted. 

We can also mimic the Minkowski solutions for excited states, although in 
the general case with explicit t dependence in all the quantities the physical 
interpretation is not immediately evident. Nevertheless, we can look for solutions 
of the form 

* x [<P,t] = Aq(t)V>i[<M (88) 
= 2N.it) J d d x^h x J d d yj^h y <P{x)G{x,mt)f {x) {y) MM 

(89) 

This solves the Schrodinger equation with the same kernel as for the vacuum 
state. The analogue of the energy gap equation is 

AE^t) J d d y^h~ y G(x,y;,t)f {x) (y) = Jgg {u i + m 2 + iR) x f (x) (x) (90) 

where 



d r, 



'dt 



In 



iVn 



(91) 



The nature of these states depends on the functions f(\)(x). These may be chosen 
to be eigenfunctions of the operator + m 2 + £R), the suffix A denoting the 
set of quantum numbers specifying the degenerate eigenfunction (the analogues 
of the momenta k in Minkowski spacetime). 

We can also write down analogues of the creation and annihilation operators 
of section 2. If we write 



a(\;t) 



J d d x^h x j d d y^h y G(x,y]t)<p(y) + 
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X) 



f(xM) 



(92) 
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and 
a f (A;t) 

then 



d d x\[^h x J d d yj^h y G(x, y; t)<j>(y) 



5 



/(A)(£) 



a(A;t)^ [<M = 



(93) 
(94) 



and 



However, in contrast to Minkowski spacetime, there is no simple particle inter- 
pretation associated with these operators for a general spacetime. 



o(A; t), a\p; t)} = 2 d d x^h x / d d yj^h y f* {x) {x)G(x,y;, t)f {p) {y) (95) 



The kernel equation (pTj ) is difficult to solve in general, so we shall specialise 
to three cases of particular interest - 'static' spacetimes, where the metric is a 
function only of the space coordinates, 'dynamic' or Bianchi type I spacetimes, 
where the metric depends only on the time coordinate, and a special class of 
conformally static metrics of sufficient generality to include all Robertson- Walker 
spacetimes. Before restricting to these special cases, however, it is of some interest 
to give a partial solution valid in general. 

This partial solution is given in terms of a solution ip(x) to the wave equation 



y/=g[D + m 2 + £R] tp(x) = 
Any G(x,y;i) satisfying the equation 



d 



d d y^h~ y G(x,y;,t)^(t,y) = -i^jgf^{t,x) 



(96) 



(97) 



is then a solution of the full kernel eq.(87). 

To see this, multiply the kernel equation by ip(t, y) and integrate with respect 
to y and then substitute (|97j). This gives 

d 



i J d d y ipfay)— [yfh x h y G{x,y\t) 



+ % J d d y^h x h y G(x,y;,t)^ [ip(t,y) 



and combining the two differential terms we find 



d 



d d y_\j-h x ^hy G(x, y\ t)tp(t, y) 



Using eq.(|97j) for the kernel again now gives 



0_ 

Ot 



x 
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which is just the wave equation satisfied by i/j(t,x). 

Later, when we have found explicit solutions for static or dynamic spacetimes, 
we see that these may be found from eq.(|97D by writing the function ip(x) as an 
appropriate Fourier transform, using eqs.( |216|) or (|220|) respectively, and using 
orthonormality and completeness relations where applicable. 



3.1 Static spacetimes 

For static spacetimes, the kernel equation is solved with a time-independent kernel 
satisfying 

/ tfV=& G(a;,2) G{z,y) = sjg^ip, + m 2 + £R) X 5 d (x,y) (98) 

The solutions can be expressed in terms of the eigenfunctions of (□j+m 2 +£i?), 
which are the Fourier transformed solutions -ip^(uj,x) of the wave equation for a 
static metric, viz. 

(□, + m 2 + £R)$ (x) (uj,x) = g 00 u 2 (X)ij (x) (u,x) (99) 

The eigenfunctions are specified by a set of discrete or continuous quantum num- 
bers A which generalise the momentum label k in flat space. The eigenvalue 
&>(A) is a function of the A. These eigenfunctions satisfy orthonormality and 
completeness relations described in appendix A. 

It is then straightforward to check that the kernel equation is solved by 



G(x,y) = J $ {x) (u,x)r w (u,y) (100) 

In passing, we also note that this expression follows by expressing ip(x) in eq. (]9"T|) 
as a Fourier transform and using the completeness relation (21S). It is also clear 
that the kernel is a real function. 

In Minkowski spacetime, we showed that the kernel was related to the inverse 
Wightman function evaluated at equal time. We can derive a similar result for a 
general static metric. Introducing the inverse of the kernel, A(x,y), defined by 

J d d z^h z A(x,z) G(z,y) = 5 d (x,y) (101) 
we find using the orthonormality and completeness relations that 

A(x,y) = J |^ -L $ w (u,x)r w (u,y) (102) 

The Wightman function is evaluated in appendix B, eq. ( |246| ) and as expected we 
confirm 

A(x,y) = 2g + (x,y)\ ET (103) 
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Again the reason for this identification is given in appendix |C[ 

The complete Schrodinger vacuum wave functional can therefore be written 

as 

= #o(*M>M (104) 

where 



fpo = exp 



(105) 

Notice how the fields 0(x) effectively project out the kernel functions ip^iu^x) 
appropriate to the boundary conditions on on the hypersurface S. We shall dis- 
cuss this point in more detail in the examples which follow. The time dependent 
phase factor is 

N (t) = e~ iEot (106) 



where E is the divergent vacuum energy, 

E = - J d d x\l-h xA [cg G(x,x) 



- I dfi(X) w(A) 5 d (0) 



(107) 
(108) 



The analysis given above of the excited states for a general spacetime may 
now be applied. For a static spacetime, all the quantities are time independent, 
so the states generated by the creation and annihilation operators are genuine 
stationary states with respect to the chosen Hamiltonian evolution. 

The wave functional for the first excited state is 

^[<f>,t] = 2iV\(t) J d d xf^h x J d d yj^h y 4>{x)G{x,y)f {x) {y) M4>] (109) 

where G(x,y) is the static kernel, Ni(t) = e -1 "^ N (t) and /m (y) is a particular 
eigenfunction of + m 2 + £R) specifying the state. The corresponding creation 
and annihilation operators simplify to 



a(A) 



at (A) 



d d x 



d d x 



V^V^ w(A)0(aO + 
^K^[° uj(\)<f>(x) - 



8 



X) 



5 



5(j)(x) 



/(A)fe) 
/(A)(z) 



and satisfy the commutation relations 

[a(A),ot(p)] = 2w(A) (2vr) d 5 d (A,p) 



(HO) 
(HI) 

(112) 
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As an elementary example, it is straightforward to check that in Minkowski 
spacetime choosing ip(\){oJ\x) — expik.x, where u(k) = \fk? + m 2 and dfx(\) = 
d d k, reproduces precisely the results of section 2. 



3.2 Dynamic spacetimes 



By a 'dynamic' spacetime, we mean one where the metric depends only on the 
'time' coordinate t. The general metric, after rescaling the t coordinate, can 
therefore be written as 

ds 2 = dt 2 - a 2 At)dx i dx j (113) 



and defines a Bianchi type I spacetime 



To simplify the kernel equation for the special case of a dynamic spacetime, 
we introduce the Fourier transform G(k,t) (see appendix A for definitions), 



G{x, t) 

The kernel equation reduces to 
. d 



d d k 
(27r) d 



ik-(x—y) 



y J G(k; t) 



(114) 



dt 



hG(k;t) =h v ^G 2 (k;t) 



-guj 2 (k; t) 



(115) 



where ui 2 (k; t) = —g^kikj + m 2 + £R. 

The solution^ may be expressed in terms of the functions 4>(t,k) introduced 
in appendix [A] which satisfy the Fourier transform of the wave equation, viz. 



We find 



G(k; t) 

and thus for the kernel itself, 



(□ - g ij kikj + m 2 + £,R)ip(t,k) = 
—i d 



-g dt 



\mjj(t,k) 



G(x, y; t) 



1 



d d k 



ik.(x-y) 



d 



lmp(t,k) 



, -h J (2n) d dt 
Notice that this is again consistent with the general expression 



(116) 
(117) 

(118) 



above. 



To show this [13[], we convert the non-linear first-order differential equation 
for G(k; t) into a linear second-order equation. Define 



F(k;t) = exp{-i dtahG(k;t)} 



(119) 



6 Thc Schrodinger formalism in fiat Robertson- Walker spacetime has been extensively studied 
in refs.Jl|, 0. 
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where a 



~S~ and let d = %. We then have 
-h at 



q2 o 

^-F = -ia^-{hG) - a 2 h 2 G 2 F 
at 1 at 



iahGF 



and using eq.( J115[ ), 



: F = -— F -a^^uj 2 F 



(120) 



(121) 



dt 2 adt 

However, this is just the equation satisfied by the functions ip(t,k), since by 
expanding D in eq. (|116| ) we see 



d 2 d + ay 

a 







(122) 



We therefore identify F(k;t) with ip(t,k) and eq.( |117|) follows. 



The function ip(t,k) is the solution of a second-order ordinary differential 
equation and is therefore a linear combination of two independent solutions with 
arbitrary coefficients. Since the kernel itself is the derivative of the logarithm of 
ip{t,k), the overall normalisation is unimportant. The kernel therefore depends 
on one arbitrary parameter. 

This is a physically important difference from the static spacetime, where 
the kernel is essentially unique. For a dynamic spacetime, the 'vacuum' states 
described by the Gaussian solution of the Schrodinger wave functional equation 
form a one-parameter family [fL3| . The selection of one of these as the vacuum 
must then be made on further physical grounds. We shall return to this point in 
the context of particular examples later. 

The Schrodinger vacuum wave functional can therefore be written as 



* [<M = N (t)MM 



(123) 



where 

and 
with 



cxp 



-h 



d d k d r 



In V>(t, k) I d d x I d d y (Kx)e^- ( ^V(?/) 



(27r) d dt 



N (t) 



1 



E (t) = -y/9o~oV-h I d d k G(k;t)5 d (0) 



-if* dtE (t) 



(124) 
(125) 

(126) 



In this case, however, Eq does not have a useful interpretation as an energy. 

The analogue of the higher excited states can be found by substituting for 
the kernel in the general formulae following eq.(B3), where fk(x) 



e 1 -- is an 
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eigenstate of + wi 2 + £R) with eigenvalue u 2 (k;t). Again, however, we em- 
phasise that the physical interpretation of these states is not as straightforward 
as for a static spacetime where they are genuine stationary states of the energy 
associated with the chosen Hamiltonian evolution. 

It is important to check the stability of the vacuum state defined by the 
time-dependent wave functional ( |124j ). Vacuum stability requires that 



V<j>%((f>,t)* (M = \N (t)\ 2 J V^oiMMM (127) 
should be time independent. From eqs. (|125|) and (|126|) we have 

|iV (t)| 2 = exp{-l Jd d xJ0- d InH 2 } (128) 

where we have used eq. (|117|) for the kernel. Obviously this will in general be time 
dependent. 

The second contribution to eq.( |127| ) is 

P</»|^o| 2 ~ [detihGnfayt))]-* 

= exp{-~| d d xj 'A^\n[hG^(k;t)]} (129) 

where is the real part of the (generally complex) kernel. From eq.( |117| ) we 



sec 



ft(fc») = fp^ (130) 

while for the imaginary part, 

e * tH v^F^' 2 < 131 > 

where W[ip, ip*} = ij) (J^*) — 4>* (JiVO is the Wronskian of the solution and its 
complex conjugate. Using eqs.( |114| ) and ( |130|) therefore gives 

J V<P \M 2 ~ exp{+i Jd d xJ^- d In |^| 2 } (132) 

up to time-independent terms. These remaining terms are seen to be time inde- 
pendent using Abel's theorem, viz. that g 00 y/—gW is a constant, where W is the 
Wronskian defined above. 

The overall time dependence therefore cancels as required, confirming the 
stability of the vacuum state in a dynamic spacetime. 
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3.3 Conformally static and Robertson- Walker spacetimes 

A straightforward generalisation of the results of sections 3.1 and 3.2 allows us 
to find the kernel for the class of (d + l)-dimensional spacetimes with metric 

ds 2 = g 00 {t)dt 2 - a 2 (t)da 2 (133) 

where da 2 is a d-dimensional static space with metric (Tij(x). Let a = deter,.,-. 
Introducing a rescaled 'conformal time' rj, the metric can be written as 

ds 2 = C{rj) drj 2 - a ij (x)dx i dx j (134) 

The spacetime is therefore conformally static, but with the conformal scale factor 
restricted to depend only on time. This class includes the important Robertson- 
Walker spacetimes, which are special cases of eq.( |133| ) where goo = 1 and da 2 is 
a 3-dimensional static space of constant curvature k = —1, 0, +1, viz. 

da 2 = d X 2 + f(x)[d0 2 + sin 2 Odcj) 2 } (135) 

where 

!sinx < x < 2vr k = +1 
X 0< X <oo k = (136) 

sinh x < x < 00 k = — 1 

For k = — 1, 0, +1 this metric describes hyperbolic (open), flat, and spherical 
(closed) spaces respectively. 

Since the spatial metric is static, general theorems guarantee the existence of 
orthonormal spatial modes y(\)(x) which play the role of the Fourier functions 
expik.x in the previous section. A full discussion of the solutions of the wave 
equation for these spacetimes is given in appendix A. 2. 

We therefore write the kernel as 

G{x,m t) = f^$ G(X; t) y {x) (x) y* w (y) (137) 



Then, using the solutions T(t, A) of the transformed wave equation (see eq. ( [225] )), 
we readily find (c.f. eqs.( |117| ) or (p7|)) 

G(A;*) = -i-^L J^lnT(t,A) (138) 

The kernel for a conformally static spacetime with a conformal factor which is a 
function of time only is therefore 

. — i f da(X) \d 
G(x,y;t)^- 1 



fgo-oa«J (2.y [Ot ^ X) \ ^ ^ 

The interpretation is identical to that for the dynamic metric case. Again, there is 
a one-parameter family of vacuum states corresponding to the freedom in selecting 
the solution T(t, A) of the wave equation. 
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4 Cosmological Model I 



To illustrate this formalism, we now consider two simple two-dimensional 'cos- 
mological models' of Robertson- Walker type. We give an exact description of the 
vacuum states, emphasising the role of boundary conditions, and discuss carefully 
the phenomenon of particle creation due to the cosmological expansion. 
The spacetime metric is chosen to be of Robertson- Walker form, 

ds 2 = dt 2 - a{t) 2 dx 2 (140) 

Rescaling the time coordinate such that dt = a(t)drj, the metric may be rewritten 
in terms of the 'conformal time' rj as 

ds 2 = C(r l )(drf -dx 2 ) (141) 

The metric is conformal to Minkowski spacetime. 



o 



A+B 






A-B 


















-infinity infinity 

n 

Figure 1: Conformal scale factor for the cosmological model with expansion. 

The first model [0J is chosen to have asymptotically static regions separated 
by a period of expansion (see Fig(|I])). The conformal scale factor is chosen to be 

C(rj) = A + B tanh(pr?) (142) 

where A,B,p are constants with A > B > 0. In the asymptotic IN (rj — > — oo) 
and OUT (rj — > oo) regions, 

C{rj) -> A ± B as rj -> ±00 (143) 

This choice allows explicit analytic solutions for the wave functionals in terms of 
hypergeometric functions. 

The metric is in the class we have called 'dynamic', so the construction of 
the vacuum wave functional is a straightforward application of the techniques 
described in section 3.2. The first step is to find the solutions of the wave equation 
in the metric specified by eq. (|142|) . 
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4.1 Wave equation 

Taking the Fourier transform, 

ip(rj,x) 

the wave equation becomes 



00 dk 



Akx 



-00 (2?r) 



[d 2 + (k 2 + Am 2 ) + 5m 2 tanh(pr/)]^(^, fc) = 



(144) 



(145) 



Notice that the only occurrence of the conformal scale factor is with the mass 
m. This is a trivial consequence of the conformal invariance of the scalar field 
equation for m = (we have taken £ = for simplicity). Obviously any non- 
Minkowskian behaviour of the vacuum must involve the mass. 
It is convenient at this point to make the definitions: 



uj 2 (k;r]) 
.2 



oj 4 . 



out 



k 2 + m 2 C{7]) 
k 2 + m 2 (A — B) 
k 2 + m 2 {A + B) 

77 (^W i ^in) 



With the substitution 

and a change of variable z 
hypergeometric equation 



[2 cosh(pr])]~~r x(v,k) 



(146) 
(147) 
(148) 

(149) 



(150) 



|[1 + tanh(p^)], the wave equation reduces to the 



z(l - z)d 2 z + 



1 



IWj. 



p 



-z(2 + 2 



iw. 



P 



a- 



p 



+ 



wt 



p z 



X(z,k) = (151) 



The properties of this equation and its solutions are discussed in [[16], |17[ . Here, 
since — 00 < 77 < 00, we are only concerned with solutions in the range < z < 1. 

Because of the singular points at z = and z — 1 in the hypergeometric 
equation we have to consider two sets of solutions to cover this range completely. 
The first set are valid in the range < z < 1, while the second set are valid for 
< z < 1. The radius of convergence of the two sets of solutions is strictly less 
than 1. 

Within the radius of convergence of the z = singular point there are two 
linearly independent solutions. We can choose these so that as z — > 0, one, ijj™, 
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behaves as a positive frequency solution while the other, tp™, becomes a negative 
frequency solution. Explicitly, 



ip™(r],k) = e- iw+v [2cosh(pr])}-— x 

2^(^,1 + ^;!-^) 



(152) 



r '^°° e - iw inV 



(153) 



e - iw+r >[2 cosher?)] ~ 



; (1 + tanh(pr/)) 



x 



^(sst.l + sstjl + isja^) 



(154) 



^^0° e +iw in rt 



(155) 



The solutions valid around the z = 1 singular point can similarly be chosen 
so that they correspond to positive and negative frequency solutions in the limit 
z — > 1. These are 



^(77, fc) = e-^pcosh^)]-— x 



2 Fi( ! 



— ,1 + — ;l + ^;l-z 

p ' p ' p ' 



(156) 



(157) 



-tw+r] 



[2 cosh(p?7)]" 



2^1 



P ' P 



— tanh(pr^)) 
■1-z) 



™out 

p 



X 



(158) 



(159) 



In the overlap region < z < 1, both sets of solutions are valid and we may 
express one in terms of the other. Specifically, we have 



^l n (77,fc) = a(k)tlj°_ ut + p(k)ijj. 



\out 



(160) 



where 



Y(\ -\- iWin )Y( iw ° ut 
a{k) = " ' ; 



r(^)r(i + ^) 



(3{k) = — '- 1 " - 



r(^)r(i-^) 



(161) 
(162) 
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The magnitudes of these two coefficients are 



\<*(k)\' 



W: 



sinh 2 (ZE^±) 



out J sinh sinh (^) 

sinh 2 (^) 



sinh f 2^ J sinh f ^ 



and satisfy 



|a(fc)| 2 -|/3(fc)p 



(163) 



(164) 



(165) 



Notice that the solution which is positive frequency in the z — > IN region is 
a mixture of the solutions of both positive and negative frequencies in the z — *■ 1 
OUT region. As we shall see, this mixing is at the heart of the 'particle creation' 
interpretation. The a and (3 coefficients above are just the Bogoliubov coefficients 
in conventional treatments J18, Q. 



4.2 Vacuum wave functional 



The vacuum wave functional is given by the standard formula ( |118| ) for a dynamic 
spacetime. That is, 



G(x,y;r)) 



00 dk 



-oo (2tt) 



Jk( x -y) G ( k . v } 



where 







ln[ai)™{r],k) + ty™{r], 



(166) 



(167) 



C{rj) drj 

To specify the vacuum beyond this one-parameter ambiguity, we need to im- 
pose a boundary condition. We choose this so that in the IN region, the wave 
functional becomes the usual positive frequency Minkowski vacuum 



(168) 



This implies we choose a = 0, b = 1 above, i.e. use the negative frequency solution 
of the wave equation in the dynamic kernel. 
Explicitly, 



G(k; V ) 







[lnC n (^)] 



C{rj) drj 

— |r ju; + + u;_ tanh(pr?) - u in [l - tanh(p?y)] 
2 F 1 (>+, 1 + >+; 



(169) 



3*1 ( 



^+,l + ^ + ;l + ^;z 



(170) 
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which in the asymptotic limit is simply 

G(k; r) — > — oo) 



(171) 



(A-B) 

as required to reproduce eq. (|168|) . The complete vacuum wave functional is there- 
fore 

%[t,v]=N (7i)M$,v] (172) 



where 



and 



I r°° dk ~ ~ 



N q (t]) = exp 



driC(rj) / dkG(k; rf)6(0) 



(173) 



(174) 



This is the full analytic expression, valid for all values of the time rj, for the 
vacuum state chosen to satisfy a positive frequency boundary condition in the IN 
region. It is stable for all times, the proof being precisely that given in section 3 
for a general dynamic metric. 

The next step is to examine the behaviour of this state in the asymptotic fu- 



ture OUT region. For this, we use eq. (160) to reexpress the kernel in terms of the 
wave equation solutions and ip™* which have simple asymptotic behaviour 
as z — > 1. This gives 



G(k; V ) 



—i 



out\ 



(175) 



As rj — > oo this function does not tend to a fixed value but rather to a limit circle 
in the complex plane. We find 



G{k]T} 



oo 



^out 



(A + B) 



7e 



-2iuioutri 



1 + 7e 



-2i(jj out ri 



(176) 



where 7 = -. The field-dependent part of the wave functional is therefore 



1 r°° dk 



00 



cxp 



2tt 



out 



7e" 



-2iu) ou tr) 



1 + -ye 



-2iuj utV 



\4>{k)\< 



(177) 



4.3 Vacuum states and particle creation 



We now discuss the interpretation of this vacuum state and in particular the 
phenomenon of particle creation. First, notice that eq. ( |1 77| ) is not the standard 



Minkowski vacuum corresponding to the asymptotic limit of the metric, viz. 

1 r°° dk 



[0] = exp 



2tt 



Uout \<j>{k)\' 



178) 
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In fact, this is the asymptotic limit of a different state ^' [4>, f]} chosen from the 
one-parameter family of vacuum solutions to the Schrodinger equation, specified 
by the kernel (c.f. eq. ( |169| )) 

G\k-,v) = -^^[lnr_ ut ( V ,k)} (179) 

Completeness of the solutions to the Schrodinger equation means that any 
state can be expressed as a superposition of one of the one-parameter family 
of vacuum solutions and the excited states built from it. In particular, we can 
express the chosen vacuum state \I/o [<^, ?7] as a linear combination of the state 
[(f), ff\ and its excited states. In the OUT region, this amounts to expressing 
the Gaussian state ipo[(j), r\ — ► oo] of eq. (|177|) in terms of the excited states of the 
Minkowski state ip$ UT [(f>\. 

In the limit 77 — > 00, these excitations are standard Minkowski particle states 
and we can characterise them by their particle number and momenta. The prob- 
lem therefore reduces to the one solved in section 2.1, i.e. expressing a general, 
time-dependent Gaussian state in Minkowski spacetime in terms of a superpo- 
sition of many-particle states. The time-dependent kernel Q is read off from 
eq. dT77|) , viz. 



n(k; T]) = u out 



(180) 



Transcribing results directly from section 2.1, we find that the expectation value 
of the number density of particles with momentum k in the state *f>o[(j),r] — > 00] 
is 

(tfofe* - ocW(k)\nlv^ 00]) = (181) 

Remarkably, despite the residual time-dependence in Q, this expectation value is 
independent of rj. Explicitly, 



1 - M ^in 

sinh 2 (g^) 
sinh f^t) sinh (^) 



183) 



This is the Schrodinger picture derivation of the traditional result on particle cre- 
ation in this expanding universe model. Notice that the particle number density 
is proportional to |/5| 2 , the Bogoliubov coefficient in the conventional treatment, 
which is only non-zero if there is a mixing between the positive and negative 
frequency solutions of the wave equation in the respective asymptotic regions. 
As expected, (Af(k)) also vanishes in the conformal limit m — >• 0. 
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It is also illuminating to evaluate the expectation value of the energy momen- 
tum tensor in these states. In particular, (^ol^ool^o) evaluated in the asymptotic 
IN and OUT regions measures the time-independent energy eigenvalue of the state 
\I/o- In the intermediate region, of course, \l/o, like all the other 'vacuum' solu- 
tions, is not an energy eigenstate. With the definition of T 00 implicit in eq.([83|) 
(recall that H = J d d x y/—gg 00 T Q0 ), we find 



<*„ M T M |* oM = -C(,)/_- G (*;,)[l - - m 

where we have normalised ^ such that ($ |^ ) — 1- 
In the IN region, we simply have 

I rod ^fc 

(^o|Too|*o)/7v = ■= 7T u ™ ( 185 ) 

This is just the usual Minkowski zero-point energy density. On the other hand, 
substituting (|176|) into eq. (|184| ) we find, after some calculation, 



2 



;i84) 



/oo Ale , \\ 

— u out UjV(k)) + -) (186) 
-oo Z7T v 2 7 

Notice that the rj dependence from the kernel G(k; rj — > oo) has again cancelled 
out. ty [<j),ri — > oo] is therefore not the minimum energy state in the Minkowski 
OUT region. Rather, eq. (|186|) describes the energy of a many-particle state with 
(J\f(k)) particles of energy u out (k) as described above. 

The Schrodinger picture therefore brings a particular clarity to the issue of 
particle creation in this model. It is evident that there is no state which satisfies 
all the properties of a standard Minkowski vacuum state. As we have repeatedly 
emphasised, there is a one-parameter family of Gaussian 'vacuum' solutions to 
the wave functional equation for this class of dynamic spacetimes. If we choose to 
consider the state ^o[4>, rj] which in the asymptotic past is the Minkowski vacuum 
state "4>q N [4>\ and follow its evolution, we find that in the asymptotic future OUT 
region it is described by the time-dependent kernel ( |176j ). At this point, there 
is really no more to be said - that is a complete description of the state. It is 
stable and any physical quantity may be evaluated by taking appropriate operator 
matrix elements. The remaining question is simply to describe what the state 
[0? ?7] looks like. As we have shown above, in the OUT region where the 
spacetime is again asymptotically Minkowskian, it is precisely a superposition 
of many-particle states built from the minimum energy Minkowski vacuum state 
4>q UJ [(f>]. In the asymptotic OUT region, therefore, it is true to say that particles 
have been created as a result of the expansion. For any intermediate value of the 
time, \&o [0, v] continues to provide a complete description of the state, although 
in this case there is no simple interpretation in terms of Minkowski particles. 
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5 Cosmological Model II 

In this section, we consider a second cosmological model of the same type H], 
this time with a conformal scale factor 

C(V) = ,f, ; (187) 
cosh (pr)) 

This is shown in Fig(|2|). This describes a universe which expands exponentially 



A A 2 




-infinity 



infinity 



Figure 2: Conformal scale factor for the recollapse model. 

from an initial singularity before eventually recollapsing. (Since C{rj) —>■ as 
i] — > ±oo, the size of the universe tends to zero in the past and future asymptotic 
limits.) Notice that although the universe is infinitely long-lived in terms of the 
conformal time rj, the range of the Robertson- Walker time coordinate t in eq.( |140| ) 
is only from t = to t = T = Air/ p. In fact, solving the relation dt 2 = C(rj)drf, 
we find a(t) = Asm pt/A, representing an oscillating universe of which we are 
considering only the first cycle. 

The analysis of this model follows closely that of section 4. Again we have two 
asymptotically flat regions in which we can compare the vacuum solutions of the 
dynamic spacetime Schrodinger equation with the equivalent Minkowski vacua. 
With some additional subtleties, the results are similar and we again find particle 
creation occurring as a result of the gravitational expansion and recontraction. 



5.1 Wave equation 

In this case, the Fourier transformed wave equation is 

[fig + k 2 + A 2 m 2 cosh" 2 (prfiltifa k ) = ( 188 ) 
Changing variables from 7] to z as before and making the substitution 

~ ih 

ip(r),k) = [2 cosh(p^)] - p xiVik) (189) 
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the wave equation reduces to the hypergeometric equation (TO, IT? 



z(l - z)dl + 



1 + 



ik 
P 



z{2 + 2 



ik, 
~P' 



d z + 



k 2 ik A 2 m 2 
p 2 p p 2 



X(z,k)=0 



(190) 

Again, there are two sets of solutions of interest, valid in the neighbourhood of 
the singular points z = and z = 1 with a radius of convergence strictly less than 
1. Choosing a basis of asymptotically positive or negative frequency solutions, 
we have 



[2cosh(pr ] )]- f 2 F 1 (f - a, 1 + f + a; 1 + f ; 



(191) 



^rZ 00 e + ikr i 



k) = = pcosh^n^r^F! (-a, l+a;l ik 



p ' I 

(192) 



and 

^°- Ut (v,k) = [2cosh(/^)]-f [l-z]-f 2 F 1 (-a,l+a;l-^;l-z 



(193) 



v l^ c e + ikr i 



r + ut (v, k) = (V-'T = [2 cosh(pr / )]-V 2 F 1 - a, 1 + f + tr; 1 + f ; 1 - z) 

(194) 



where 

o = ~\ ± ^1 + ^ (195) 

Notice that these solutions are independent of the mass m in the asymptotic 
limit. The mass decouples when the size of the universe shrinks beyond the scale 
of vrC 1 . 

In the intermediate region < z < 1, we can express these solutions in terms 
of each other. In particular, we will need the relation 

k) = a{k)^°_ ut + f3{k)^j° + ut (196) 



30 



where 



a{k) 



(3(h) 



r(i + f)r(f; 



r(f - <x)r(i + f + a) 
r(i + f)r(-f) 



The magnitudes of these coefficients are 

sin 2 (7ro") cosh 2 (7r&;) + cos 2 (7rcr) smh 2 (nk) 



(197) 
(198) 



\a(k)\ 2 

\m\ 2 



sinh (irk) 



SUl 7TO" 



and satisfy 



sinh 2 (7r/c) 

\a{k)\ 2 - \(3{k)? 



(199) 
(200) 

(201) 



5.2 Vacuum Wave Functional and Particle Creation 

The construction of the vacuum wave functional is identical to section 4. Choos- 
ing boundary conditions such that the vacuum wave functional \l/o[0, v] becomes 
the standard Minkowski vacuum in the asymptotic IN region, viz. 



k \$(k)\ l 



(202) 



fixes the one-parameter freedom of vacuum states. This choice determines the 
kernel, 

i d 



C(rj) drj 



[In fc)] 



(203) 



Explicitly, 

G(k; V ) 



if 1 
■ , , \ -iktanh(pri) + -[1 -tanh 2 (p?7)] x 



:f -*)(! + £ + cr) 2Fl (l + ^-a,2 



ik 
P 



a:2 + ^:z 



2^1 (f-^l 



P 



204) 



Rewriting eq. ( j203| ) in terms of the wave equation solutions ij)™* valid in the 
OUT region, we have 



G(k; V ) 



a{ip°_ ut ) + 



(205) 
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Then, taking the asymptotic future limit r] — ■> oo, we find the following expression 
for the dependent part of the vacuum wave functional: 



oo 



exp 



dk 

W) 



k 



1 — 7(1 — z) p 
1 + 7(1 - zft 



(206) 



where 7 



This is to be compared with the equivalent Minkowski vacuum in 



the OUT region, 



= exp 



dk 
-00 (2tt) 



fc l?(*)| S 



(207) 



which is of course identical to ipQ N in this case. The same argument as in section 
4.3 now shows that the asymptotic vacuum state ( |206| ) is a linear superposition 
of many-particle states built from the Minkowski vacuum ipQ UT . The expectation 
value of the particle number density is 

(*o[0,»7->oo]|JV(A:)|*o[5,^->oo]) = 



(208) 





1 


2 


1 - 




7 


2 



Sin 7T(J 



sinh (jrk) 



in agreement with the canonical approach ||. Once again, therefore, we see 
that the initial Minkowski vacuum has evolved in the asymptotic future to a 
state described by a Gaussian wave functional with a time-dependent kernel, 
which may be interpreted as a many-particle state. We conclude that despite 
the recontraction following the initial exponential expansion, particles have been 
created in this cosmological model as a result of the time-dependence of the 
metric. 

A remarkable feature of this model, however, is that for certain values of 
the parameters, there is no particle creation. From eq.( 208 ) we see that (Af(k)) 
vanishes when a is an integer n, i.e. 



Am 
P 



n(n + 1) 



(209) 



Equivalently, there is no particle creation when the cycle period of the Robertson- 
Walker metric and the mass are related by 



T = — \/n(n + 1) 7r 
m 



(210) 
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6 Outlook 



These simple models already illustrate some of the power of the Schrodinger 
picture techniques to describe vacuum states in spacetimes with time-dependent 
metrics. As we have seen in the discussion of cosmological particle creation, 
there is a clear advantage in such spacetimes of representing the vacuum state by 
a Gaussian wave functional, entirely characterised by a simple kernel function, 
rather than the usual Fock space description. This approach also makes clear 
that it is only in the very special case of certain static spacetimes that 'vacuum' 
states exist which share all the properties associated with the usual Minkowski 
vacuum. 

In the companion paper, we extend this investigation of the Schrodinger pic- 
ture to simple model spacetimes with boundaries. In this case, even when the 
spacetime is static, the Unruh effect sharpens the question of defining a vacuum 
state by requiring us to identify which, if any, class of observers (modelled in 
the Schrodinger picture by the choice of foliationP]) perceives the chosen state to 
resemble a Minkowski vacuum. 

A complete understanding of all these issues is a necessary precursor to de- 
scribing vacuum states and particle creation in Schwarzschild-Kruskal black hole 
spacetimes, which comprise both static and dynamic regions separated by bound- 
aries which are event horizons. 
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A Wave Equation 



The wave equation for a massive scalar field is 

(□ + m 2 + £R)ip(x) = (211) 
where the Laplacian operator is 

□ = -Ljjd^y/^gdv) (212) 

and g = det(g^ u ). The wave equation can be rewritten as 

(□ + nii + m 2 + Ci2)V>(aO = (213) 

where 

□< = -j=d t {g^^d 3 ) (214) 
□o = -LdoGTv^o) (215) 



-g 

in coordinates where the metric components g^ vanish. 
A.l Static and dynamic spacetimes 

If the metric is independent of a coordinate, then there exists a Killing vector 
associated with the corresponding translation symmetry. We use this to intro- 
duce convenient Fourier transforms for the special cases of static and dynamic 
spacetimes. 

In the static case, because of the translational invariance of the wave equation 
with respect to t, the solution may be written as 

— -e— ty(",flO (216) 

-co [ZTT ) 

where 

(□< - g°°u 2 + m 2 + (R) x) = (217) 

This equation has degenerate solutions labelled by a discrete or continuous set of 
quantum numbers (A), which generalise the momentum k in flat space. 

An extensive discussion of the general properties of such solutions in static 
spacetimes has been given by Fulling [19|]. They satisfy the orthonormality 
and completeness relations 

d d x 



^/^V^W w >^Vto(w,aO = *(A,p) (218) 



34 



and 

/ f^F $w( u >v) = sf&P&v) (219) 

where the measure //(A) and delta function S(X, p) are such that 

dp(X) 5(X,p) /(A) = /(p) 



In the dynamic case, because of the translational invariance of the wave equa- 
tion with respect to x l , the solution may be written as 

= /(0? e "^^) ( 22 °) 

where k J , the measure is defined in terms of the covariant vector 

components as d d k = dk\ . . . dk d and ip(t,k) satisfies 

(□ - g ij kikj + m 2 + £R)ijj(t, k) = (221) 

A. 2 Conformally static and Robertson- Walker spacetimes 

Now consider spacetimes with a conformally static metric of the form (|133|) . The 
Laplacian n Q becomes 

□o = ^a-%(^a d do) (222) 
and all dependence on the spatial coordinates drops out. The Laplacian { is 



□ i = -^di^y/ddj 



-1 

o?\fa 

a~ 2 U d (223) 



where n d is the Laplacian on the d- dimensional static space with metric <Jij. The 
wave equation takes the form 

[a 2 U - n d + a 2 (m 2 + £R)]ip(x) = (224) 

The solutions factorise and we find 

^(t,x) = y w (x)T(t,X) (225) 

where y(\)(x) are eigenfunctions of O d with eigenvalues K 2 (X), i.e. 

O d y w (x) = -K 2 (X)y w (x) (226) 
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and T(t, A) is a solution of the transformed wave equation 

[a 2 a + K 2 {\) + a 2 {m 2 + £R)]T{t, A) = 



(227) 



All the analysis above for static spacetimes applies to eq.( |226| ). General theorems 
assure the existence of a complete, orthonormal set of solutions y( X )(x) satisfying 



d d x 



(228) 
(229) 



Robertson- Walker spacetimes are a special case of this class of conformally 
static metrics. The solutions of the wave equation for k = —1, 0, +1 are therefore 



pV^^^^d) = s(x, P ) 



where 3^a are eigenf unctions of the 3-dimensional Laplacian D 3 , 

a 3 y x ( x ,e^) = -K 2 (\)y x (x,o,<p) 

and T(t, A) satisfies 

[a^dt^dt) + K 2 (X) + a 2 (m 2 + £R)]T(t, A) = 



(230) 
(231) 
(232) 



Here, K 2 (X) = k 2 — K, where (A) = (k, j, m) and the eigenf unctions y^ X ) (x) , |20|, ^TJ 
are explicitly 



3W 



X 



(2n)m ( kj \x)Y- 



(2n)2 J±J lH (k X )Y™(e,. 



K = -1 
K = 
K = +1 



(233) 



The Y m are spherical harmonics and the function IT ) is defined by 



nfc } (x) = [^k 2 (k 2 + l)...(k 2 + j 2 )]-f sinh' X 



I (P +1 cos{kx) 
\G?(cosh xy +1 



(234) 



The function Uj^ is obtained from IlL by replacing k with — ik and x with —%X- 
Alternative forms of this decomposition are given in 0, [22| . 



The orthonormality and completeness relations are given by eqs.( |228| ) and 
( p9|) , viz. 

f^^y w (x)y* {p] (x) = 5(X,p) 



-j^y(X)(x)y w (y) = - 



(235) 
(236) 
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where ^fa = / 2 (x) sin#. The integration measures for the different types of 
Robertson- Walker metric are 



J d/j,(X) 



Jo dk Ys'jLo J2m=-j K — 



(237) 



which also specifies the form and range of the quantum numbers k, j, and m. 

For the spatially flat k = metrics, it is of course frequently more convenient 
to choose the eigenfunctions in Cartesian form 

y W {x) = exp(^) (238) 

where (A) = (hi, k 2 , k 3 ) and K 2 (\) = k 2 . The measure is just / d/i(X) = j d 3 k. 

B Green Functions 

The Green function is given by the equation 

+ m 2 + £R)Q{x, y) = (a; - y) (239) 

As in the previous section, if the equation is translationally invariant with respect 
to a particular coordinate then the Green function may be expressed as a Fourier 
transform. 

In the static case, 

Q(x, y) = r -gL e-^-'^g^; x, y) (240) 

where 

v^(n, - g 00 uj 2 + m 2 + £R)g(u;x,y) = -5 d (x - y) (241) 
For dynamic spacetimes, 

Q(x, y) = J-0j- d e ik -^G{k x°, y°) (242) 

where 

v^(D - 9 ij hkj +m 2 + ZR)Q(k; x°, y°) = -5(x° - y°) (243) 

We cannot say more about the dynamic Green function in general. However, in 
the static case, we can give a general solution. Using the orthonormality relations 
of solutions to the wave equation, we find 

dw f dfi(\) e-^o-w) $ w ( u , x) r w (v,y) 



G ^ y) = Uw)JW (^#) (244) 
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plus solutions of the homogeneous equation. The 00 integral is along the real axis, 
but there are poles at w — ±z/. The choice of contour (ie prescription) selects 
one of a variety of different Green functions [I]. [2[] . 

For example, the contours giving the Wightman function are shown in Fig(|3|). 
Evaluating the integral gives 

g + (x,y) = (0| (pix) cp(y) |0> (245) 

= ll^h^^^^^^ (246) 

As the Wightman function is just the difference between two Green functions, it 
satisfies the homogeneous wave equation. 
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Figure 3: The omega contour integral for the Wightman functions. 



C Expectation values and two-point functions 

The vacuum expectation value of an operator 0(ip, it) (a function of the field ip 
and conjugate momentum 7r) in the Schrodinger picture is 

<0| 0{<p, tt) |0) = jv<p% 0(0, -i— ) $ (247) 

The action of the momentum is simply 

= i J tfyyfh^^G^mt) *o (248) 

To evaluate the action of the field, it is convenient to modify the wave functional 
by formally introducing a 'source', viz. 

= #0 exp { l - J d d x^K J{x) <f>{x)} (249) 
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Omitting field independent terms for clarity, the vacuum amplitude is then 



(0|0)j ~ j V(J)e~ f dd - $ ^yV^y^-^^-'K^^K) + 1 d d xs/-h x J{x)<t>(x) (250) 

~ exp Au(x,y;t)J(y)\ (251) 

Notice that the functional integral involves just the real part G$ of the kernel. 
Also note that Asr is defined here as the inverse of G^ (not the real part of the 
inverse kernel). Expectation values of products of fields can then be found by 
repeated differentiation with respect to the source, e.g. 



(0\<p(x)<p(y)\0) = V<f>^*j<f>(x)cf>(y)^j 



15 15 



(0|0>. 



j=o 



V^K SJ(x) JZh y 5J(y) 
-A n (x,y;t) (252) 



This is the equal-time Wightman function, confirming the identifications (f20|) 
and ( |103| ) in the text. 

Other useful two-point functions are 

(0\n(x)n(y) |0) = -^Jhjiy{ G^(x,yj, t) 

+Jd d u Jd d v/hJ^Gz(x,u- t)A dt (u,v;t)GQ(v,y^,t)}(253) 

and 

(0\[<p(x),ir(y)}\0) = i5\x-y) (254) 
(0\{<f(x),7r{y)}\0) = - J d d u^/hJ^GQ{x,u;t)A m {u,y : ,t) (255) 
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